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Abstract —  In  this  paper,  we  present  new  iterative  algorithms  to  find  a  root  of  the  given  nonlinear 
transcendental  equations.  In  the  proposed  algorithms,  we  use  nonlinear  Taylor ’s  polynomial  interpolation  and 
a  modified  error  correction  term  with  a  fixed-point  concept.  We  also  investigated  for  possible  extension  of  the 
higher  order  iterative  algorithms  in  single  variable  to  higher  dimension.  Several  numerical  examples  are 
presented  to  illustrate  the  proposed  algorithms. 
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I.  INTRODUCTION 

Solving  nonlinear  transcendental  equations  is  one  of  the 
important  research  areas  in  numerical  analysis,  and  finding 
the  exact  solutions  of  nonlinear  transcendental  equations  is 
a  difficult  task.  In  several  occasions,  it  may  not  be  simple 
to  get  exact  solutions.  Even  solutions  of  equations  exist; 
they  may  not  be  real  rational.  In  this  case,  we  need  to  find 
its  rational/decimal  approximation.  Hence,  numerical 
methods  are  helpful  to  find  approximate  solutions. 
Numerical  iteration  methods  for  solving  nonlinear 
transcendental  or  algebraic  equations  are  useful  in  applied 
sciences  and  engineering  problem.  A  nonlinear  equation  in 
n  variables  does  not  have  the  linear 
form  f  =  anxn+  an_xxn_x  +  •  •  •  +  axxx ,  where  an , an_x ,...,ax 

are  constants,  and  also  it  does  not  satisfy  the  conditions 
such  as  f(ab)  -  f(a)f(b)  or  f{a  +  b)  -  f(a)  +  f(b) .  For 

example,  the  transcendental  equation  z  =  xex  whose 
solution  is  given  by  the  Lambert  W  function,  W(z),  defined 

as  W(z)ew{z)  =z ,  which  is  a  fixed  point  equation  of  the 

form  p(x)  =  x .  This  model  equation  has  several 

applications  in  science;  Physics,  Biology,  Material  science, 
Quantum  statistics,  differential  equations,  signal 
processing,  in  the  theory  of  zeta  functions  etc.  [1].  Very 
especial  examples  of  algebraic  nonlinear  equations  are  the 

r\ 

polynomials  Pn  (x)  =  a0  +  axx  +  a2x  H - h  anxn  .  Solving 

science  and  engineering  problems  is,  directly  or  indirectly, 
a  root  finding.  In  general,  solving  transcendental  equations 


is  a  more  difficult  when  compared  with  algebraic 
polynomials,  for  example,  the  polynomial  solving  problem 
can  be  transformable  to  linear  eigenvalue  problem  solving. 
The  complexity  of  the  polynomial  problem  depends  on 
their  degrees  and  number  of  terms. 

A  root-finding  problem  is  to  solve  for  a  root  r  also  called  a 
zero  or  a  solution  of  the  equation  /(x)  =  0 ,  such  that 

f(r)-  0.  The  root  finding  methods  are  available  in  the 

different  forms,  for  example,  the  methods  can  be  in  the 
form  of  direct/symbolic,  graphical  or  numerical  iterative. 
Solving  transcendental  equations  for  multiple  roots  is  more 
difficult  than  solving  for  a  simple  root  using  iterative 
methods.  The  available  iterative  methods/  algorithms  are 
derived  from  interpolations  (linear,  quadratic,  cubic), 
perturbation  method,  variational  technique,  fixed -point 
methods,  series  estimations  and  so  many  others.  The 
Taylor’s  interpolation  technique  has  also  many 
applications;  solving  differential  equations,  function 
estimations  and  computing,  developing  and  analysis  of 
numerical  schemes,  etc.  Taylor’s  method  is  a  foundation 
for  calculus  of  finite  differences. 

Given  any  equation  can  be  rewritten  as  an  equivalent 
equation  g(x)  =  p(x)  with  the  same  root  r.  The  iteration 

method  to  estimate  the  root  r  is  then  g(xn+x)  =  p(xn) .  The 

iterative  method  converges  when  I  p  '(x0 )  /  g  \x0 )  l<  1 ,  with 

an  appropriate  initial  guess  xo  [2].  In  the  literature,  there 
exist  several  root  solving  methods  such  as  Bisection 
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method,  Secant  method,  Regula  Falsi,  Newton’s  method 
and  its  variants  methods  to  accelerate  convergence,  namely 
Chebyshev’s  method,  Halley’s  method,  Super  Halley’s 
method  etc.  [3-12].  These  methods  may  be  classified  as 
fixed  (single)  point  methods  or  methods  based  on  choice  of 
intervals.  Nevertheless,  choices  of  initial  guesses,  interval 
selections,  and  existence  of  derivatives,  acceleration 
convergence,  efficiency  and  possible  extensions  to  two- 
dimensions  or  three-dimensions  are  some  common 
drawbacks  in  connection  with  algorithmic  complexities  [6- 
24].  In  the  present  paper,  we  apply  Taylor’s  series 
approximation  of  a  nonlinear  equation  /(x)  combined 

with  fixed-point  concept  and  some  modified  error 
correction  term  to  derive  new  iterative  methods.  The 
concept  of  function  estimation  is  applied.  The  function 
estimation  concept  is  an  extension,  improvement  and  a 
modification  of  works  in  [21],  see  also  [22,  23].  Extensions 
of  the  proposed  algorithms  to  higher  dimensions  are  also 
considered  to  solve  the  given  nonlinear  systems.  Various 
interpolation  techniques  and  root-finding  algorithms  are 
available  in  [25-33]. 

The  paper  is  organized  as  follows:  Section  2  presents  the 
proposed  algorithms;  in  Section  3,  we  discuss  the 
convergence  analysis  of  the  proposed  algorithms;  in 
Section  4,  we  present  numerical  example  to  illustrate  the 
algorithms  also  presented  comparisons  with  other  existing 
methods.  Finally,  Section  5  presents  conclusion  and 
recommendations. 

II.  PROPOSED  NEW  ITERATIVE  ALGORITHMS 

In  this  section,  we  derive  fixed-point  iterative  functions  of 
the  form  x  =  'F(x)  from  which  iterative  models 

xn+l  =  0(xn )  are  obtained. 


Consider  the  Taylor’s  polynomial  interpolation  of  a 
nonlinear  function  /(x)  differentiable  about  a  number  xo 

as: 

f{x)  =  f<x„)  +  f(x„)  +  f  -(%)  <x  -*«  >2 

n\ 

with  x  is  a  root  r  of  /(x) ,  very  small  error  |Ax-|  -h  , 
where  r-xG  -  h-  A x{ .  Now 

f(x0  +  Axi)  =  f(x0)  +  A xj  '(x0 )  +  0.5  Ax:,2/  "(x0 )  +  •  •  •  ( 1 ) 

Iff  =  f(x,y),  then  the  Taylor’s  series  expansion  in  two 
dimension  is 

f(x  +  h,  y  +  k)  =  /(x,  y)  +  {hfx  +  kfy ) 

1  (2) 
+  -(h2fxx+2hkfxy+k2fyy)  +  - 

and  if  /  =  f(x,y,z )  in  3-dimension,  then 
f(x  +  h,  y  +  k,  Z+ 1)  =  f(x,  y,  z)  +  ( hfx  +  kfy  +  //. ) 

+  \ (h2L  +  2 hkfxy  +  2 hlfxz  +  2 klfyz  •  k2fyy  I /  /.. )  +  •  •  • 

where  h  -  Ax,  k  =  Ay,  l  =  Az  are  small  step  sizes  in  x,  y ,  z 

respectively.  The  linear  estimate  from  (1), 

f(xo )  +  A xf  \xo )  «  0  gives  the  Newton’s  method  for 

simple  roots  [18]. 

Xk+ 1  =xk-  - ,  k  =  0,1,..., n.  (3) 

/  (xk) 

From  a  modified  Newton’s  method  of  order  two  for 
multiple  roots  [16],  one  has  a  second  order  convergent 
method  for  finding  simple  roots  [22] . 


Definition  1:  A  point  p  is  said  to  be  a  fixed-point  of  /(x) 
if  and  only  if  f(p)  =  p  . 

The  stipulation  of  the  fixed-point  concept  plays  very 
important  role  in  root  finding  algorithms.  For  every  root  r 
of  /(x) ,  there  exists  an  auxiliary  equation  p(x)  such  that 

f(r)  =  r-  p(r)  =  0  .  That  is  r  must  be  a  fixed-point 

of  p(x)  .  Solving  for  a  root  of  /(x)  =  0  by  finding  a  fixed- 

point  of  x  =  p(x)  guarantees  an  iterative  method 

xn+l=P(xn),  n  =  0,1,2,....  In  literature,  one  can  find 

several  methods  that  having  fixed-point  forms.  For 
examples,  Newton’s  method,  Halley’s  method, 
Chebyshev’s  method  (see  equations  (3),  (10),  (11))  etc. 


xk+ 1  ~  xk 


/(**) 

u 


k  =  0,1,..., ?z, 


(4) 


where  U  -  f  \xk )  -  f(xk )  and  h  -  - 


f 


u  +  f 


is  error/ 


correction  part  of  (3);  whereas  h'  = - is  the  error  term 

U 

of  (4).  These  error  terms  play  important  role  to  investigate 
the  proposed  methods. 

(i)  A  combination  of  (3)  and  (4)  gives  higher  order 
method  as  a  fixed-point  formula 


O(x)  =  x 


A/(x)  |  Bf(x) 
U  +  f  U 
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(ii)  If  we  choose  the  quadratic  model  of  (1),  then 
f  (xa  +  h)  =  f  (x  )  +  hf  \xQ)  +  h  f  }X°^  *  0  .  We  obtain 


h  = 


-2  /(x0) 


2  f\x0)  +  hf"(x0) 


and 


h  =  - 


/(x0)  +  Q.5/z2/"(x0) 

/'(•*, 0) 

/ 


(6) 


Using  the  error  term  h  =  —  from  (4)  in  (6),  one  can  get 


two  iterative  methods  as  follows 

2f(xk)\U(xk)\ 


xk+ 1  =  xk  ~ 


2  f\xk)[U(xk)]-f(xk)f\xk) 


(7) 


xk+ 1  =  xk  ~ 


f(xk) 


-0.5 


_  f(xk)  f"(xk) 
U(xk)  +  f(xk)  ^  f\xk)[U(xk)f 

We  have,  at  a  simple  root,  f  '  ~U  .  Hence 


(8) 


xk+ i  =xk~ 


f(xk ) 


U(xk)  +  f(xk) 


-0.5 


f(xk)2f\xk) 


[U  (xk  )]- 


(9) 


/ 


Inserting  the  usual  error  term  h  =  —  —  from  (3)  in  (6),  one 
has  Halley’s  method  (10)  and  Chebyshev’s  method  (11) 


xk+ 1  =xk~ 


xk+ 1  =  xk  ~ 


fXXf'W) 

(10) 

[/  '(x*:  )]2  —  1  /  2/ (xj. )/  "(xk ) ' 

/(**) 
/  '(x, ) 

l[f(xk)]2f'Xxk) 

2  [/'(x,)]3 

(11) 

Also  we  have  the  extension  of  Newton  method  (12)  and 
Euler  method  (14),  See  [5,  21]. 


(p(x)  =  x  - 


/’W+V[/'U)]2-2/(x)/"(x) 


/ "(*)  /  "M 

2H(x) 

1  +  ^jl-2p(x) 
H{x)f\x) 
f  \x)  ’ 

f(x) 


(12) 


(^(x)  =  x  - 


2C*)  = 


W(x)  = 


t/+/ 


(13) 


(iii)  Using  U  =  /  '(x)  -  /(x) ,  from  the  linear  estimation  of 
/(x)  in  (1),  one  obtains 


fi(x, )  =  xk  - 


/(x,)3+/(x,)U2 

U(f(xk)2+f\xk)U) 


(14) 


For  =  0,1,...,^, 


**+i  =xk~ 


f{xk)  f(xk)  +  f2(xk) 


U(xk)  +  f(xk)  f'2(xk)  +  f2(xk) 


=  xk- 


f(xk)  +  (u(xk )  +  f(xk ))' 


f 


X 


(M(Xj.)  +  /(X|t))2 

/(**)  ,  f(xk) 


(15) 


-(■ 


■y 


\u(xk)  +  f(xk)  u(xk)  +  f(xk)  j 

(iv)  Suppose  consider  the  cubic  model  of  /(x)  as 

f(x0+h)  =  f(x0)  +  hf(x0) 

+  \h2f\x0)  +  \h2f"Xx0)*  0. 

2  6 


(16) 


(16)  =>  hf '+  \  h3f  ’’’  =  h(f'+l/6h2f "’)  =  -(/  + 1  /  2/?2/  ’’) 

6 

xt  u  ,  6/+3/J2/”  T 

Now  we  have  h  - - - - .  Inserting  a  correction 


6/’+/iz/ 


f 

term  h-  —  ,  we  have 
U 


h  =  -  3 


2jT2-4/'/2+2/3+/2/" 
6/'3-12/  '2  /  +6/'/2  +  f 2  f 


fff 


This  gives  an  iteration  function  using  (i),  (ii),  (iii). 


yr(x)  =  x  — 3  2^2-4/'/2+2/3+/2/" 
6/ ,3-12/ '2  /  +  6/'/2  +/2/ 


(18) 


Replacing  the  / by 


/'/" 

/ 


,  we  get  the  method 


yr(x)  =  x-3  2#'2-4/'/2+2/3  +  /2/" 
6/  ,3-12/  '2  /  +  6f'f2+jf'f" 


(19) 


(v)  Now,  we  derive  the  new  iterative  functions  using 
Taylor  expansion  and  iterative  methods  are  given  (22), 
(23),  (24),  (25)  and  (26). 

Consider 

f(x0)  +  hf  \xa )  +  \h2f  \x„ )  +  \h2f  "\x0 )  *  0,  (20) 

2  o 


which  gives 

hf'(x0)  +  U2f"(x0)  = 

hence 


/(xo)  +  ~/  6ft3/ "'(x0) 
6 


and 
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h  — 


6  f(x0)  +  h3f'\x0) 
6f\x0)  +  3hf\x0) 

f  , 


If  we  substitute  h-  —  in  (21),  we  get 

U 


xk+ 1 =xk~2 


U[f(xk)-l/6U3f  "\xk ) 


2 f\xk)U  -  f{xk)f\xk) 
xk+i=xk  -2 (/'(x*)-/(xj) 

(  s,..  ^  A3 


X 


/CO- -7 


/(•O 


V 


f(xk)-f(xk) 


f  ’"O ) 


k>  J 


/  [2/  \xk )(/  \xk  )  -  /  (xt ))  -  /(X, )/  "(X, )] 
Utilizing  / ^  J  in  (23),  we  have 

**+i  =xk  -2 (f\xk)-f(xk)) 


x 


/(-O-- 


/2(0 


f  ’(X,  )/  "(X,  ) 


6  (/  \xk )-  f(xk ))' 

/  [2/  \xk )(/  ’(x, )  -  f(xk ))  -  /(xt )/  "(xt )] 


/ 


ti2 


And  if  f = - ,  we  obtain 

/ 

*t+,  =xt-2(/'(xt)-/(xt)) 


X 


/(**)“ 


0(0 


/  "O ): 


6(/'(xt)-/(xjy 
/  [2/  '(xt )(/  ’(Xt )  -  /(xt ))  -  /(x, )/  "(X, )] 
One  may  also  obtain 

„  UVf\xk) 

Xu ,  1  -  Xu-  2 - - . 

2  f\xk)U-f(xk)f\xk) 
where  U  =f\xk)- f(xk )  and 


V  =  f(xk)~ 
6 


f(xk) 

u 


(21) 


(22) 


(23) 


(24) 


(25) 


(26) 


III.  CONVERGENCE  ANALYSIS 

In  this  section,  we  present  the  convergence  analysis  of  the 
proposed  methods.  The  basic  concepts  and  definitions  can 
be  found  in  [2,  4,  5,  18]  for  further  details. 

Definition  2:  Let  ( xn )  be  the  sequence  of  approximate 

roots  obtained  using  the  iterative  method  xn+l  -  p(xn )  for 

exact  root  r.  The  convergence  is  said  to  he  reached  if 
lim  p(xn )  =  r  ,  or  equivalently  lim  ( en )  =  0.  The  sequence 

of  errors  is  denoted  by  en  -  r- p(xn ) . 

Theorem  1:  (Fixed-Point  Iteration  Theorem)  Let  f(x)  =  0 
be  written  in  the  form  x  =  p(x) .  Assume  that  p(x)  has 
the  following  properties: 

•  For  all  v  in  [ a ,  b],  p(x)  is  in  [ a ,  b],  i.e.,  p(x) 
takes  every  value  between  a  and  b. 

•  p  \x)  exists  on  (a,  h)  with  the  property  that  there 
exists  a  positive  constant  k  <  1  such  that 

p\x)  I  <  k  for  all  v  in  (a,  b). 

Then 

(i)  there  is  a  unique  fixed  point  p  of  p(x)  in  [ a ,  b\. 

(ii)  for  any  x0  in  [ a ,  b],  the  sequence  xn+l  -  p(xn) ,  n 

=  0,  1,  2,...  converges  to  the  fixed  the  point  p,  i.e.,  to 
the  root  r-p  of  f(x)  -  0 . 

Theorem  2:  (Order  of  Convergence)  Assume  that  (j){x ) 
has  sufficiently  many  derivatives  at  a  root  r  of  f  (x)  .  The 
order  of  any  one-point  iteration  function  (j){x )  is  a 
positive  integer  p ,  more  especially  (f){x )  has  order  p  if  and 

only  if  (j){r)  =  r  and  =  0  for  0  <  j  <  p, 


Note:  (Nonlinear  Systems  in  Two-dimensional  Space) 
The  extension  of  the  methods  in  (7)  and  (8)  can  be 
expressed  as  follows.  The  equation  (7)  can  be  written  in 

two-dimensional  space  as  xk+l  =xk  -(A)[1-IF]_1,  where 
N  =  xk-  <p(xk ).  and  IV  =  0.5  ff "(/  V1  [/ ,/T' ;  and  the 
equation  (8)  as  xi+1  =  (xk -^(x))[l- 0.5#^ "((/'- fY'fl 

Also  (9)  is  expressed  as  xk+l  =xk -(A)[1-1V(/')3] . 
Where  cp(x)  is  Newton’s  iteration  function  to  solve 
f(x)  =  0  in  one  variable.  The  Newton’s  method  to  solve  a 
system  of  equations  F(2O  =  0  is  given  by 

AXk  = -J~l(Xk)F(XK).  Where  J  is  the  Jacobean  of  F 
and  J~l  is  its  inverse  when  exists. 


<!>(p\r)*  0. 

Now  we  can  prove  the  order  of  convergence  of  algorithm 
in  (7)  as  follows.  Consider  the  iterative  method 

xk+\  =xk -(A0[1-M]_1,  or 

xk+1=xk  -{xk  -<p(xkW-0.5jrf'-l[UYlTl,  (27) 

where  N  =  xk-cp{xk),M  =0.5  jf"  f'~l\UYl  ■  We  can 
write  iteration  equation  (27)  as  </>{x)  —  x  —  (x  —  (p(x))w , 

where  w(v)  =  [1  -M[U]~l]~{  and  x  =  (p(x) ,  Newton’s 
iteration  function. 

Let  r  be  a  simple  root  of  f(x)  =  0.  We  have 
q>(r)  =  r ,  cp(x)  =  x  and  </>(x)  =  x.  And 
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but  0.  Differentiating  ,  where 

(j){x)  =  x-(x-  <p(x))w ,  we  find  that 
(f V(r )  =  (jf\r )  =  0  but  ^  0  .  Therefore,  p  >  3  . 

Conversely,  if  p  =  3,  then  we  can  show  that 
{ j)\r )  =  —  0  but  ^  0  .  Hence,  (7)  or  (27) 

is  third  order  convergent  method.  Similarly,  we  can  do  for 
the  other  methods. 

IV.  NUMERICAL  EXAMPLE 

Consider  the  following  equations  to  test  the 
efficiency  of  the  proposed  methods. 

(1)  fl(x)  =  5/2x3 -5/2x-5/2  =  0,  with  x0  =  1,2 

and  root  r  ~  1.324718  in  (1,  2), 

(2)  f2 (x)  =  6x - 2cos x  —  2  =  0,  x0  =  0,  1 ,  2, 

r~  0.607102  in  (0,  1). 

(3 ) f3(x)  =  COSX  +  X3  —  e*  =0,  x0=-2.5,  -1,  -0.5, 
r  ~  -0.649565  in  (-1,0). 

(4) /4(x)  =  4x6 -4x-4  =  0  ,  x0=l,2, 


r  ~  1.134724  in  (1,2) 

(5)  /5  (x)  =  log10  (x)  -  2x  +  2 ,  x0  =  0.5,  1.5,  2.2, 
r  =1.000000  in  [1,2)  [5]. 

The  comparisons  of  the  results  obtained  by  Newton 
method  (NM),  Chebyshev’s  method  (CM),  Halley’s 
method  (HM)  with  the  proposed  methods  in  equations  (4), 
(7),  (8),  (23),  and  (24)  are  given  in  Table  1.  The 
implementation  of  the  proposed  algorithms  in  C++  is  done, 
and  the  number  of  iterations  taken  to  converge  to  a  root  r 
with  tolerance  Tol  =  10'7  is  recorded.  The  “Functions  (f)” 
in  table  refers  to  the  number  of  functional  evaluations  up 
to  derivatives  and  “Nar”  denotes  average  number  of 
iterations.  The  efficiency  indices  can  be  calculated  as 

e  =  fllp  ,  where  p  is  order  of  convergence,  for  example,  if 

/=  p  =  2,  then  e  -  1.4142,  also  e  -  1.4142  for  p  -f-  4. 
Hence,  the  higher  order  doesn’t  give  guarantee  for  the 
faster  convergent  or  better  efficient. 


Table  1:  Numerical  Results  of  Convergence  Analysis 


f,x 0 

Nm 

(4) 

CM 

(7) 

(8) 

HM 

(23) 

(24) 

(25) 

fv  1>  2 

5,5 

3,  3 

4,3 

2,2 

2,4 

3,3 

2,3 

3,  3 

3,2 

f2  :0’  l’2 

3,3,  3 

5,5,- 

3,3,3 

2,3,3 

2,  3,3 

2,2,3 

3,2,3 

3,2,3 

3,2,3 

f3  :-2.5,-l,-.5 

7,5,4 

6,5,4 

4,3,4 

4,3,3 

4,  3,3 

4,3,3 

4,  3,3 

4,  3,3 

4,3,3 

/4X2 

5,7 

3,5 

3,4 

2,2 

2,2 

3,4 

2,3 

2,3 

3,2 

f5  5, 1.5, 2.2 

3,3,  3 

4,5,6 

3,2,3 

3,2,3 

3,2,2 

3,2,3 

2,  2,3 

2,  2,3 

2,2,3 

Functions  (/) 

2 

2 

3 

3 

3 

3 

4 

3 

3 

Nar 

4 

4 

3 

3 

3 

3 

3 

3 

3 
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for/](x):  Green  (Newton  method).  Red  (7).  Blue  (23) 

Fig.l:  Graphs  of  three  fixed  point  methods 


V.  DISCUSSIONS 

In  Table  1,  the  average  number  of  iterations  elapsed  to 
converge  to  the  root  r  and  numbers  of  function  evaluations 
up  to  derivatives  have  been  determined.  The  graphs  of  the 
fixed-point  functions  for  some  methods  are  also  shown. 
The  use  of  the  error/correction  term  results  in  fast 
convergent  methods,  but  may  not  affect  the  number  of 
functional  evaluations  up  to  derivatives.  Moreover,  it  can 
be  seen  from  the  graphs  that  the  new  methods  converge  at 
single  fixed-point  (root  r)  even  for  larger  intervals 
containing  the  root  r  when  compared  with  existing 
method.  It  is  also  found  that  higher  order  does  not  mean 
better  efficient.  Efficiency  e  depends  on  number  of 
function  evaluations  /  and  order  p.  In  addition  to  new 
methods,  the  work  has  improved,  in  terms  of  proofs, 
graphs  and  results,  by  few  authors.  In  this  study,  we  also 
reviewed  several  applications  of  root  finding  algorithms  in 
science  and  engineering.  For  example,  (1)  to  find 
intersection  of  curves,  (2)  computing  the  critical  or 
stability  points  of  functions,  (3)  to  determine  extreme 
values  of  optimization  problems,  (4)  for  solving  matrix 
eigenvalue  problems  and  applicable  in  other  designing 
problems. 

VI.  CONCLUSIONS  AND  RECOMMENDATIONS 

In  the  paper,  we  applied  Taylor’s  nonlinear  polynomial 
interpolation  and  a  correction  term  with  fixed -point 
concept  to  obtain  new  iterative  methods  for  estimating 
simple  roots  of  nonlinear  transcendental  or  algebraic 
equations.  The  correction  technique  applied  to  both 


quadratic  and  cubic  models  does  affect  convergence.  We 
have  shown  possible  extensions  for  solving  higher 
dimensional  nonlinear  systems.  In  our  future  work,  we 
shall  present  further  analyses  of  these  algorithms  and  other 
higher  order  iterative  algorithms  with  applications. 
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